Efficient electron spin manipulation in a quantum well by an in-plane electric field 
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Electron spins in a semiconductor quantum well couple to an electric field via spin-orbit interac- 
tion. We show that the standard spin-orbit coupling mechanisms can provide extraordinary efficient 
electron spin manipulation by an in-plane ac electric field. 
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Efficient manipulation of electron spins by an exter- 
nal ac field is one of the central problems of semicon- 
ductor spintronicsi, quantum computing and informa- 
tion processing^ The original proposals of spin manipu- 
lation were based on using a time dependent magnetic 
field B(t). However, there is a growing understand- 
ing of the advantages of spin manipulation by a time- 
dependent electric field E(t) that couples to the elec- 
tron spin through different mechanisms of spin-orbit (SO) 
interaction^ Recently Kato et al. successfully manipu- 
lated electron spins in a parabolic AlGaAs quantum well 
(QW) by a gigahertz bias applied to a single gate; this 
produced a field E{t) perpendicular to the well— The 
structure was engineered to make the electron Lande ten- 
sor g small (|(/| < 0.1), highly anisotropic, and position- 
dependent, and to achieve in this way a rather strong SO 
coupling originating from the position-dependence of the 
Zeeman energy. The present authors have shown£ that 
for large g-factors typical of narrow-gap A3B5 semicon- 
ductors such a field E(t) can provide efficient electrical 
operation via the standard mechanisms of SO coupling 
(Dresselhaus- and Rashba^) but only for relatively wide 
QWs. Indeed, the coupling of 2D electron spins to a per- 
pendicular field E(t) develops due to a deviation from 
the strict two-dimensional (2D) limit and is proportional 
to w/X (for w/X <C 1), where w and A are the confine- 
ment and magnetic lengths, respectively. Using a tilted 
magnetic field B is critical for the mechanisms of Refs.0 
andyj 

In this paper we show that using an in-plane electric 
field E(t) results in a significant increase in the coupling 
of this field to electron spins. Moreover, the coupling ex- 
ists in the strict 2D limit (and usually increases with de- 
creasing w) and for arbitrary orientations of _B, including 
a perpendicular-to-plane B. However, using a tilted field 
B allows distinguishing the contributions coming from 
the competing mechanisms of SO coupling. Of course, 
producing an in-plane field E(t) needs designing a proper 
coupling to the radio frequency or microwave source; this 
problem can be solved. 

We consider a 2D electron gas in a A3B5 crystal con- 
fined in a (0,0,1) QW subject to a tilted magnetic field B. 
The 2D kinetic momentum of electrons k = —i'V+eA/hc 
depends only on the normal component of B, B z = 
BcosO, hence A(r) = (B z /2)(—y,x,0). The opera- 
tors k x and k y obey the usual commutation relations 



[k x ,k y ]- = —i/X 2 , where A = (ch/eBz) 1 / 2 ; we assume 
that B z > 0. The total Hamiltonian H = H + H z + 

H so + H int includes the orbital term Hq = h k /2m, 
where m is the electron effective mass, the Zeeman term 
Hz = g/J-sicr • B)/2, /j,b being the Bohr magneton, and 
the coupling Hi nt (t) = e(r-E(t)) to an in-plane ac electric 
field E(t). The SO interaction H so = Hr> + Hr includes 
the Dresselhaus and Rashba terms that in the principal 
crystal axes are 

Hd = a D (a x k x - (Tyky), H R — a R (a x k y - a y k x ). (1) 

It is convenient to change from the set of operators 
{x, y, —id x , —id y ) to the kinetic momenta (k x , k y ) and the 
coordinates of the center of the orbit, (xq, yo)~ The latter 
ones commute with (k x ,k y ) and with Hq and obey the 
commutation relations [xq, yo\- = iX 2 . As usual, instead 
of (k x ,ky) the Bose-operators a — X(k x — ik y )/^/2 and 
a + = X(k x + iky)j\pl can be used. After these transfor- 
mations, we come to the convenient expressions for the 
coordinates that appear in the operator H[ nt (t) 

x = xq — iX(a + — a)/ V2, y = yo — X(a + + a)/ V2. (2) 

Without the SO interaction, the energy spectrum of 
2D electrons is described by two sets of Landau levels 
E a (n) = hcu c (9)(n+l/2) + huj s a for two spin projections 
on the B direction, a — ±1/2. Here lu c (Q) = eB z /hc = 
loj_cos9 and uj s — g^iBB/fi are the cyclotron and spin 
flip frequencies, respectively, and n > 0. The energy 
spectrum of a 2D system with a single term in H so , ei- 
ther Hr> or Hr, can be found analytically but only for a 
perpendicular field B. The problem with both terms in 
H so cannot be solved analytically. We consider in what 
follows the limit of a strong magnetic field (that is of prin- 
cipal physical interest and allows solving the problem for 
an arbitrary direction of B) and assume that the frequen- 
cies 10 c and u s are large compared with the SO coupling, 
H so <C huj Cl hu; s . Then the SO term can be eliminated, 
in the first order in H so , by a canonical transformation 
exp(T)i£ The operator T is non-diagonal in the orbital 
quantum number n and its matrix elements are 

(n',cr'\f\n,a) = {ri ,a'\H so \n,a) l[E a ,{n') - E a {n)]. (3) 

After the transformation, the time independent part of 
H conserves the spin projection on the magnetic field. 
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The operator r = (x, y) is diagonal in the spin indices, 
and Hi nt (t) drives spin- flip transitions due to the level 
mixing produced by H so . After the T-transformation, r 
acquires an anomalous part r so = [T, r] _ that drives spin 
transitions. The matrix elements of r so diagonal in n are 

(nf| ?"so |n |) = -{oj c (n | \[H so ,r] + \n I) 

+ uj s {n]\[H so ,rUn[)}/n{u 2 c -L0 2 s )- (4) 

the subscripts + and — designate anticommutators and 
commutators, respectively. Because H so includes only 
the operators (a, a + ), the operators (io,J/o) appearing 



in Eq. (|2J) make no contribution to the diagonal in n 
matrix elements of Eq. Q). These matrix elements can 
be calculated explicitly for Hp and Hr. 

In the quantum limit, when only the lowest Lan- 
dau level n = is populated, for a magnetic field 
B = _B(sin 9 cos tp, sin 9 sin tp, cos 9) and an in-plane elec- 
tric field E(t) polarized at an angle tp to the x axis, 
E(t) — E(t)(costp,sintp,0), the matrix elements l D R = 
(0 t l^soCosV' + y so smtp\0 1)d,r 01 the electric dipole 
spin resonance (EDSR) are 



II _ a D 
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[(uj c cos 9 — u) s ) sva(tp + tp) — i(uj c — lu s cos 9) cos(tp + tp)], 



[(lo c cos 9 + w s ) cos(<^ — tp) + i(w c + w s cos0) sin(iy9 — 



(5) 
(6) 



We assume that cos 9 > 0. When cos 9 < 0, one should 
also change the sign of u> c , hence, matrix elements are 
even with respect to the reflections in the (x, y) plane. 
For n ^ 0, the anticommutator in Eq. Q acquires a fac- 
tor (2n + 1) while the commutator does not depend on 
n. Therefore, Eqs. J5j and JHJ can be generalized by sub- 
stituting lu c — > (2n + l)w c , lu s — * oj s in the numerators. 

Let us first discuss the angular dependence of these 
matrix elements. The Hamiltonian of the Rashba in- 
teraction Hr is an invariant of the group Coov of the 
continuous rotations about the z axis. Therefore, l R is 
isotropic with respect to the joint rotations of E and B 
and depends only on the difference ((f — tp). However, 
this dependence is rather strong and the sign of the az- 
imuthal anisotropy depends on 9, Fig. la. For the Dres- 
selhaus interaction, Z-p is anisotropic, Figs, lb and lc. 
When ip. — > tp + tt/2 and tp — > tp + n/2, l^ D changes sign. 
As a result, the intensity of EDSR that is proportional to 
II 2 

l D shows a four-fold symmetry with respect to the joint 
rotations of B and E about the z axis, Fig. lc. Gen- 
erally, the intensity is proportional to (l D + l R ) 2 , hence, 
both contributions interfere and the intensity shows only 
the two-fold symmetry in accordance with the symmetry 
group C% v of the Hamiltonian H so . When Hp and Hr 
are of a comparable magnitude, as was found for GaAs 
QWsjSiiSiAiiiSi the effect of the interference is very strong, 
see Fig. Id. Therefore, the azimuthal dependence of the 
intensity of EDSR is a powerful tool for measuring the 
ratio of the coupling constants, aR/ap. 

It is seen from Eqs. (jSJ and JSJ that the EDSR driven 
by an in-plane field E should be seen for any direction of 
B including perpendicular to the QW plane. However, 
using a tilted field B provides special advantages. First, 




FIG. 1: Angular dependence of the EDSR intensity 
l(0,tp,tp) (x \ll + l ll D \ 2 for a (0,0,1) QW calculated from 
Eqs. JSJ and (|SJ for: (a) - an = and tp = 7r/4; (b) - 
cir — and tp = 7r/4; (c) - oir = and tp = tp; and 
(d) - an = Qfo and tp = tp; arbitrary units. The ratio 
of the frequencies uj a /u± — —0.17 as in InAs. The pole in 
I(6,tp,tp) at u>c(fl) ~ w s was cut off by adding an imaginary 
part iT = 0.Ulo± to u c (6). For large 9 values, in the vicinity 
of the pole and close to 7r/2, the figure provides only qualita- 
tive patterns of 1(6, tp, tp) because of its high sensitivity to the 
phenomenological line width F and the restrictions imposed 
by the condition hio c S> H BO - 



it allows distinguishing the contributions of the different 
SO coupling mechanisms. Second, because of the poles 
in the denominators at lu c (9) w u> s , it can allow strong in- 
crease in the EDSR intensity. A similar resonance in the 
EDSR intensity when the frequency of an electric dipole 
transition approaches ui s is known in the spectroscopy of 
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local centers liiiii Near the pole, where the denominator 
of Eq. (0} vanishes, perturbation theory in the interac- 
tion H so fails and the sharpness of the resonance is cut 
by the avoided level crossing and by the level width. 

Let us estimate now the magnitudes of the matrix ele- 
ments of Eqs. and © for in-plane EDSR, E ± z, and 
compare them to the matrix elements for the different 
mechanisms of spin- flip transitions. By the order of mag- 
nitude, l D ~ ajj/hujc and w an/huj c for B \\ z. With 
typical values of an ss czr « 10 -9 eV cm, m « 0.05mo, 
and B w 1 T, we get l D « l\ « 10~ 5 cm. A similar 
length for the electron paramagnetic resonance (EPR) is 
Iepr « |5|Ac/4 « 10- 10 cm for | 5 | « 10; A c = h/m c 
being the Compton length. Hence, ifj,/|j 3> ?epR: an d 
EDSR strongly dominates over EPR. 

In the perpendicular geometry, E || z, the characteris- 
tic length for the Rashba interaction is ~ anui s /hu}Q, 
where ujq « Ti/mw 2 is the confinement frequency. The 
factor oj s appears as a result of the electron confinement 
in the electric field direction; its presence is required by 
the Kramers theorem^ The ratio of the characteristic 
lengths ~ lo c u s /uj 2 w {mg/2mo){w/ A) 4 . There- 

fore, under strong confinement conditions, EDSR 
in the in-plane geometry is much stronger than in the 
perpendicular geometry because in the latter geometry it 
develops due to the deviation from the 2D regime^ 

For the Dresselhaus interaction, it is instructive to 
compare the EDSR intensity in the in-plane geometry to 
its intensity in 3Dii4 Usually the 3D Dresselhaus length 
is about l D D « 5/hio c \ 2 . Here 8 is the constant of the 
bulk inversion asymmetry of A3B5 compounds and is re- 
lated to old as an ~ 5 /w 2 pi^ The ratio of the character- 
istic lengths is about l D /l D D ~ (A/u>) 2 . There is a spe- 
cial 3D geometry, E || B, when l D D is especially large, 
l D D w 5/h~ui s \ 2 ; indeed, u> s /oj c is usually numerically 
small. However, in a similar 2D geometry, E || B _L z, 



the length l" D is also large, l' D ss ao/hujs, as follows 
from Eq. JSJ for cos 6* = and uj c — 0. Hence, again 
l D /l D D ~ (X/w) 2 . Therefore, under the conditions of a 
strong confinement, w <C A, in-plane EDSR in a QW is 
stronger than in 3D. Absence of the potential confine- 
ment in the direction of E and a strong confinement in 
z direction are highly advantageous for strong EDSR. 

The most important quantity characterizing the spin 
operation efficiency by a resonant electric field E{t) is 

the Rabi frequency n R = eEl/ti. With I « 4, 4 ~ 10 ~ 5 
cm as estimated above, we find that £Ir « 10 s -1 in an 
electric field as small as only about E w 0.6V/cm. This 
estimate shows that the electron spin manipulation by an 
in-plane electric field should be highly efficient. 

Electron heating by the electric field and spin relax- 
ation can hamper electrical operation of electron spins. 
However, electron heating is suppressed by a strong mag- 
netic field because, according to the Drude formula, it 
decreases with B as (u> c t p )~ 2 , t p being the momentum 
relaxation time. The spin relaxation is universal, i.e., it 
does not depend on the specific mechanism of spin op- 
eration. The above theory has been developed for the 
2D limit but the qualitative conclusions related to the 
high efficiency of in-plane operation are valid also when 

In conclusion, we have demonstrated an extraordinary 
high efficiency for the electron spin operation in quantum 
wells by an in-plane ac electric field. Spin coupling to 
an in-plane ac electric field is many orders of magnitude 
stronger than to an ac magnetic field, and is also much 
stronger that the coupling to a perpendicular-to-well ac 
electric field. 
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